Important marine areas for endangered African penguins before and after the crucial stage of moulting

The population of the Endangered African penguin Spheniscus demersus has decreased by > 65% in the last 20 years. A major driver of this decrease has been the reduced availability of their principal prey, sardine Sardinops sagax and anchovy Engraulis encrasicolus. To date, conservation efforts to improve prey availability have focused on spatial management strategies to reduce resource competition with purse-seine fisheries during the breeding season. However, penguins also undergo an annual catastrophic moult when they are unable to feed for several weeks. Before moulting they must accumulate sufficient energy stores to survive this critical life-history stage. Using GPS tracking data collected between 2012 and 2019, we identify important foraging areas for pre- and post-moult African penguins at three of their major colonies in South Africa: Dassen Island and Stony Point (Western Cape) and Bird Island (Eastern Cape). The foraging ranges of pre- and post-moult adult African penguins (c. 600 km from colony) was far greater than that previously observed for breeding penguins (c. 50 km from colony) and varied considerably between sites, years and pre- and post-moult stages. Despite their more extensive range during the non-breeding season, waters within 20 and 50 km of their breeding colonies were used intensively and represent important foraging areas to pre- and post-moult penguins. Furthermore, penguins in the Western Cape travelled significantly further than those in the Eastern Cape which is likely a reflection of the poor prey availability along the west coast of South Africa. Our findings identify important marine areas for pre- and post-moult African penguins and support for the expansion of fisheries-related spatio-temporal management strategies to help conserve African penguins outside the breeding season.

Marine ecosystems are under severe pressure from anthropogenic activities 1 , with negative impacts being observed throughout the trophic web 2 . Since 1950, the size of seabirds' populations are estimated to have decreased by c. 70% 3 and consequently, seabirds are one of the most threatened groups of birds globally 4,5 . They face threats on land and at sea, such as resource depletion due to climate change or competition with fisheries, introduced predators and diseases, breeding habitat loss, and mortality linked to fisheries bycatch [5][6][7][8] . No-take zones and Marine Protected Areas (MPAs) have been proposed as the most appropriate strategies to mitigate threats associated with prey availability 2,7,9 . However, in-depth knowledge of how a species uses its environment during different life-history stages is required at appropriate spatio-temporal scales to inform conservation strategies [10][11][12] .
While breeding, the foraging range of seabirds is constrained around a central location, as they must return regularly to incubate eggs and care for dependent young. These foraging constraints are relaxed outside of the breeding season, when seabirds often extend their distribution and target foraging habitats outside of their breeding foraging range 13,14 . The larger distribution of non-breeding seabirds and their increased amount of time

Results
Over the study period, 81 pre-moult and 19 post-moult African penguin foraging trips lasting > 20 days were recorded from Bird Island (19 and 6 individuals, respectively; 2012-2015) and Dassen Island (50 and 13 individuals, respectively; 2012-2019) and Stony Point (12 pre-moult individuals; 2018-2019 Table 1). Distributions were significantly different between stages, colonies and years, with penguins from Dassen Island travelling the farthest during both the pre-and post-moult stages compared to penguins from Bird Island and Stony Point.
Pre-and post-moult distribution of African penguins. During pre-moult trips, penguins from both Dassen Island and Stony Point mostly travelled south-eastward, with core ranges within the vicinity of their colonies and east of Cape Agulhas. In contrast, post-moult birds from Dassen Island mostly travelled north of St Helena Bay (Fig. 1). Pre-and post-moult birds from Bird Island remained close to the colony, but the post-moult core range was larger than the pre-moult range (Fig. 1). Permutation tests revealed that pre-and post-moult distributions were significantly different at both the core (54% UD) and distributional ranges (90% UD) of the Bird Island and Dassen Island penguins (Table 2, Fig. 1). In addition, there was significant interannual variability in the distribution of Dassen Island and Bird Island pre-moult penguins as shown by the permutation test ( Fig. 2; Table 2).
Pre-moult penguins dispersing from Stony Point were found to return to Stony Point to perform their annual moult. In contrast, 50% of the pre-moult penguins dispersing from Dassen Island, for which moulting could be determine, moulted at Dassen Island and 50% moulted at Stony Point. However, the location of moult for penguins could only be determined for penguins which were tracked with GPS-GSMs tracking devices (see "Methods" for more detail).
The pre-moult path metrics were significantly different between colonies (perANOVA: p < 0.01 for all path metrics), with penguins from Dassen Island and Stony Point travelling significantly further than penguins from Bird Island ( Fig. 3; Supplementary Table S1). Path metrics were also significantly different between pre-and postmoult penguins from Dassen Island (perANOVAs: path length: p = 0.01; maximum distance: p = 0.03; individual core areas: p = 0.28) and Bird Island (perANOVAs: path length: p = 0.98; maximum distance: p = 0.68; core areas: p = 0.02; Fig. 3). Pre-moult penguins from Dassen Island travelled significantly longer distances (path length) and significantly further from their colonies (maximum distance) than post-moult penguins from Dassen Island ( Fig. 3; Supplementary Table S1). In contrast, for penguins from Bird Island, path lengths and maximum distances were similar between pre-and post-moult trips, whereas individual core areas were significantly greater for post-moulters compared to pre-moulters ( Fig. 3; Supplementary Table S1).
Identification of marine Important Bird Areas. All three pre-moult IBAs and the Dassen Island postmoult IBA (Fig. 4) had representativeness scores > 85% (Table 3). Existing no-take zones for purse-seine fisheries, i.e. within MPAs, had little overlap with these IBAs, except for the Bird Island pre-moult IBA where there was a 32% overlap ( Table 3). The IBAs overlapped substantially with the proposed 20 km no-take zones around Dassen, Robben, St Croix and Bird islands, but not Stony Point (Table 3). Despite this variable overlap between IBAs and 20 km no-take zones, birds spent a substantial proportion of time within 20 km (16-67%) and 50 km (30-89%) of their colonies during both pre-and post-moult stages (Table 3).

Discussion
Successful conservation management of seabirds' marine habitat should incorporate all life history stages of the seabirds. We show that although pre-and post-moult adult African penguins have much larger distributions compared to during the breeding season and disparate distributions during the pre-and post-moult stages, they use waters within 20 and 50 km of their colonies intensively. This suggests that the 20 km no-take zones around colonies, proposed to alleviate competition for resources between breeding African penguins and fisheries, will also benefit the penguins outside of the breeding season. Our research highlights the need for stage-specific conservation and management of mobile species. The pre-and post-moult adult African penguins travelled up to 600 km away from their colonies, greatly exceeding the maximum foraging range of breeding penguins(< 50 km) 41,56 . Similarly to the initial dispersal of juvenile African penguins 43 , the pre-and post-moult adults rarely exceeded the 200 m depth contour but did not range as far as the juveniles. Sherley et al. 43 showed that juvenile African penguins from Dassen Island went northward of their colony to approximately 20°S whereas all pre-and post-moult Dassen Island adults in this study remained below 29°S. The maximum extent travelled along the west coast was similar for pre-moult adults and juvenile penguins dispersing from Stony Point (< 30°S). However, whereas the core areas for the Stony Point juvenile penguins was along the west coast 43 , the core area of the Stony Point pre-moult penguins was along the south coast, east of Cape Agulhas. Together with the findings of Sherley et al. 43 , we have shown that areas north of St Helena Bay on the west coast, in the vicinity of the Stony Point colony and east of Cape Agulhas are important to penguins outside of the breeding season (i.e. the pre-and post-moulting and juvenile stages) and should be adequately reflected in Marine Spatial Planning initiatives within South Africa.
Travelling to distant foraging areas during the pre-and post-moult adults 19,30,[63][64][65][66][67][68][69][70][71][72][73] and juveniles [74][75][76] stages is known for other penguin species. Due to the aggregation of high densities of breeding seabirds foraging from a central location, prey availability is thought to decline outward from these central locations as the breeding season progresses, likely leading to increased competition for available prey 77,78 . These large migrations to distant foraging habitats by the African penguin outside of the breeding season may be a relic of this mechanism, wherein the penguins disperse to distant and likely more profitable foraging areas once they are no longer responsible for dependant young (i.e. they are no longer constrained around a central location).
The pre-moult distribution of African penguins departing from Dassen Island and Bird Island showed significant inter-annual variability. Plastic foraging distributions during the pre-and post-moult stages has     [79][80][81][82][83][84] . The marine environment is highly dynamic with numerous biophysical factor determining the spatiotemporal distribution of prey 85 . Behavioural plasticity, such as plastic distributions in seabirds, may be an adaptive trait to locate ephemeral prey 86 . This is likely the case for pre-moult African penguins as plasticity in foraging behaviour as a response to variable prey distribution has been demonstrated for this species during the breeding season 87,88 .
Pre-moult African penguins in the Western Cape (Dassen Island and Stony Point) travelled significantly further than those in the Eastern Cape (Bird Island). Most of the Western Cape pre-moult penguins travelled south-eastward along the South African coastline and had core distributional ranges east of Cape Agulhas. The south-eastward movement of these birds mimic the relatively recent south-eastward shift in distribution of their main prey: anchovy and sardine, again highlighting the possible plastic foraging distribution of these birds during this stage. Due to compounding effects of environmental changes and intensive fishing pressure along the west coast, the distribution of the two forage fish species have shifted south-eastward, east of Cape Agulhas, with sharp decreases in sardine abundances along the west coast [89][90][91] . The significantly greater distances travelled and areas covered by Dassen Island pre-moult penguins compared to those from the others colonies indicates  www.nature.com/scientificreports/ that this may be as a result of shifting prey distributions and as a consequence they expend the most energy to fatten up during this crucial life history stage, compared to birds at other colonies. This may compromise their ability to accumulate sufficient body reserves for their annual moult, which may influence post-moult survival or recovery, and subsequently, pre-breeding condition and breeding success. The Dassen Island population has shown one of the fastest rates of decrease since 1999 47 and the survival of African penguins from west coast colonies has been linked to regional estimates of spawner biomass of sardine in the preceding non-breeding season 23,39,44 . Our results provide support for these findings and highlight the importance of prey availability to these birds during the pre-and post-moult life-history stages. African penguins from Bird and Dassen islands used different areas during their post-moult trips compared to pre-moult trips. However, the distribution of post-moult birds from both colonies overlapped with areas of upwelling and associated cool sea surface temperatures and high chlorophyll-a concentrations. The core ranges of post-moult penguins from Bird Island extended north of Port Alfred, along the eastern edge of the Agulhas shelf where regular upwelling is associated with relatively high prey biomass 92 . Post-moult penguins from Dassen Island mostly travelled to an area north of St Helena Bay, which prior to the south-eastward shift in the distribution of anchovy and sardine, was an important spawning area for sardine 89,93,94 . During summer, when the post-moult penguins were tracked, there is an increase in wind-driven upwelling at both of the sites which promotes phyto-and zooplankton productivity [95][96][97][98][99] , both important prey items of sardine and anchovy 100,101 . African penguins are attracted by olfactory cues, i.e. dimethyl sulphide that is released by phytoplankton 102 . Higher concentrations of such olfactory cues during summer, both within the St Helena Bay and Port Alfred areas, may explain the attraction of penguins to these areas during the post-moult stage (November/December). Alternatively, the penguins may target other cues that attract them to the cold sea surface temperatures and high chlorophyll-a concentrations typical of upwelling in these regions and also attract Cape gannets Morus capensis 42 and juvenile African penguins 43 . Previous tracking studies of pre-moult 30,63,64,70,71 and post-moult 19,[66][67][68][69][70]72,73 penguins have shown that many species have strong associations with environmental variables (e.g. chlorophyll-a concentrations) or physical features (e.g., oceanographic fronts) where their prey are concentrated.
Post-moult body condition may also have some bearing on the large disparity in foraging distributions of pre-and post-moult African penguins from Dassen Island. Due to the energy demands of their catastrophic moult, African penguins are in poor body condition at the onset of the post-moult period 32 . This may impede their ability to travel to distant foraging grounds, such as the Agulhas shelf habitat visited during their pre-moult stage, especially if there is limited prey available close to their moulting sites to improve body condition prior to departure to distant foraging areas. It may be that the penguins drift more or less passively to the area north of St Helena Bay within the jet current that flows northward along the continental shelf 98 . However, Sherley et al. 43 found that juvenile African penguins, moving in the same direction and to similar areas as the Dassen Island post-moult adults, were swimming actively to reach this area. Further fine-scale investigation into drivers of adult penguin movements outside of the breeding season is needed.
Marine Important Bird Areas for pre-moult African penguins from three of their six largest colonies 47 , and post-moult African penguins from Dassen Island were identified. All four of these newly identified areas were highly representative of their colonies' populations (as determined by their representativeness scores, see "Methods" for more details). A small proportion of these IBAs overlap with existing no-take zones for purse-seine fisheries. However, large proportions of the Dassen Island and Bird Island IBAs overlapped with the proposed 20 km no-take zones around four major colonies 103    www.nature.com/scientificreports/ 20 km and 50 km of an African penguin's colony are used intensively by pre-and post-moult penguins. As the breeding season of the African penguin is protracted and asynchronous with moulting and breeding individuals being present in the colonies throughout most of the year 108-110 , improved foraging conditions close to their colonies could have appreciable benefits for African penguins year-round. Given the endangered status of the African penguin, linked to their rapidly declining population 47 , these areas should receive high conservation priority. However, the much more extensive distribution of pre-and post-moult and juvenile 95 African penguins compared to that of breeding penguins highlights the need for life history stage-specific spatial management and conservation strategies [10][11][12] .

Conclusions
The congruent population declines of the African penguin, Cape gannet and Cape cormorant Phalacrocorax capensis in South Africa have been attributed to insufficient availability of high quality prey or increased foraging effort to reach profitable foraging areas, resulting in lower reproductive success and survival [40][41][42][43][44] . A key strategy to conserve these endangered seabirds should be to increase opportunities for the birds to access these highly mobile prey while they are available in important foraging areas during all life history stages. No-take zones around major colonies have been proposed as a strategy to improve African penguin breeding success 24,52,53 .
Our findings suggest that such fishing closures may also improve prey availability to penguins outside of the breeding season, potentially improving their survival and breeding propensity. However, future research should investigate oceanographic drivers of the penguins' pre-and post-moult distribution as this will help disentangle the differential foraging distributions between important life-history stages and elucidate the relative influences of natural versus anthropogenic drivers of prey availability for these birds. Marine spatial planning, in the form of MPAs or no-take zones, has been identified as one of the most appropriate tools to manage prey depletion by fisheries [5][6][7][8] and these measures should be expanded to incorporate important foraging areas of African penguins outside of their breeding season, such as the pre-and post-moult stages and the initial dispersal of juvenile birds.

Methods
Data collection. The lifecycle of the African penguin is relatively asynchronous and protracted compared to other penguin species as breeding and moulting penguins are often present in the colonies year-round 109,110 . However, in general, breeding peaks during the winter months (February-September and January-July at colonies to the west and east of Cape Agulhas, respectively [108][109][110] ) with the majority of birds moulting during early summer from September to January [108][109][110] . Using re-sightings of flipper bands, the pre-and post-moult stages of African penguins have been estimated to be c. 35 and c. 42 days in length 109,111 , with birds either returning to their breeding colonies or colonies closer to their pre-moult foraging grounds to moult 112 Table 1). From August to October, loggers were deployed on pre-moult breeding adults identified as those with late stage chick(s) ready to fledge (i.e. chicks with few to no downy feathers) 109 . During November, loggers were deployed on adults in the last stages of their feather moult (i.e. few to no old feathers). Loggers were attached to the feathers on the dorsal midline of the bird's lower back using tesa ® tape (Beiersdorf AG, Germany) and secured with cable ties and cyanoacrylate glue (Loctite 401 ® ). Due to the lower battery life of the PTT and Catlog GPS devices used between 2012 and 2015 (Table 1), these devices were scheduled to only record locations at night, between 18:00 and 24:00 GMT and 21:00 and 22:00 GMT, respectively, to maximise the devices' ability to log a position while the penguins rested on the surface during night-time hours. The greater battery life of the solar-powered Pathtrack GSM-GPS devices, used between 2017 and 2019, allowed for locations of the penguins to be recorded every hour. Pre-and post-moult distribution of African penguins. There was considerable variability in the number of days (defined as the full 24 h of each calendar date) individuals were tracked (range 2-78 days). Number of days tracked were defined from the first location recorded at sea to either the last location recorded at sea before evidence of moulting was noted (> 48 h on land) or the last transmitted location at sea. To ensure that tracks were representative of pre-and post-moult foraging trips, only tracks > 20 days were analysed. This is because it was found that for penguins which moulting could be determined, all pre-moult foraging trips were greater than 20 days (range 21-57 days; Table 1) Penguins roosted on land regularly (both at the deployment colony and at other colonies), therefore, locations on land were removed and tracks were split into trips between these land-based events. Tracking data were filtered for erroneous fixes based on transit speeds of greater than 12.4 km·h −1113 using the R package argosfilter 114  www.nature.com/scientificreports/ For each tracked bird, the following path metrics were calculated: (1) the maximum distance travelled away from the deployment colony, (2) the path length from the deployment colony to the maximum distance using the sum of great circle distances between consecutive locations and (3) the area covered by their core range (further referred to as 'individual core area'). These individual core areas were estimated as utilisation distributions UD115 using the R package adehabitatHR 116 with a smoothing factor (h) of 7 km following Dias et al. 58 . The isopleth that demarcated each individual core area was estimated based on optimal isopleth value selection (OIVS), following Vander Wal and Rogers 117 . The OIVS method uses the exponential relationship between the proportion of home range area used by an individual and the isopleth volume to identify thresholds (slope = 1) delineating areas of maximum use. The OIVS was applied to each individual and the mean optimal isopleth value of all individuals (54% UD) was taken to represent the core range of the birds and used in all individual core area estimations.
Due to different sampling regimes between 2012 and 2015 (locations only recorded at night) and 2017 and 2019 (1-h intervals), all tracks were down-sampled and linearly interpolated to have one location per day prior to the calculation of path metrics. Using data collected between 2017 and 2019, we investigated the impact of down-sampling the tracks by testing for significant differences in path metrics calculated with the original (24 locations per day) and down sampled data (1 location per day) using a permutational analysis of variance test (perANOVA, 5000 permutations). Maximum distance (p = 0.87) and individual core area (p = 0.42) were found to be similar between the sampling regimes. However, path length (p < 0.01) was significantly longer when sampled more frequently. Despite this bias, we retained path length because it was highly correlated between the two sampling regimes (Spearman's correlation coefficient = 0.93). In addition, the influence of the type of tracking device on path metrics and number of days tracked was investigated using a perANOVA. Due to the relatively low sample of different devices used for different stages and colonies ( Supplementary Fig. S1 online), comparisons were limited to Dassen Island's pre-moult stage. Maximum distance (p = 0.43), individual core area (p = 0.14) and path length (p = 0.26) were found to be similar between device types however number of days tracked was significantly greater for birds equipped with PTTs compared to birds equipped with GPS-GSMs and Catlogs. However, the majority of the pre-moult trips tracked with GPS-GSMs were complete even though the number of tracked days by these devices were lower than that of PTTs and Catlogs (Table 1, trips were labelled as complete if there was evidence of the bird remaining on land for > 48 h before data stopped being transmitted and could thus be assumed to be moulting). Data from all tracking devices were therefore retained. Whether the birds were moulting could not be determined for those tracked by PTTs and Catlogs as these devices were only set to record at night which prevented us from distinguishing whether birds were roosting on land overnight or moulting. Path metrics appeared to be greater for post-moult Dassen Island birds that were tracked with PTTs compared to those tracked with GPS-GSMs ( Supplementary Fig. S1 online). These data were still pooled and used in further analyses as between year comparison were not made.
Permutational analysis of variance tests (perANOVA, 5000 permutations) were also used to assess path metrics differences among colonies during the pre-moult stage. Due to the lack of post-moult data from Stony Point, differences were not investigated among colonies during the post-moult stage. Differences between preand post-moult path metrics were tested separately for Dassen Island and Bird Island. Dunn's post-hoc tests with Bonferroni corrections for multiple comparisons were applied to all significant permutation test results (package: dunn.test) 118 .
Individual penguin data were then pooled per stage and year (hereafter referred to as data group) to investigate per-stage and annual differences in core (54% UD) and distributional (90% UD) ranges. Comparisons were limited to data groups that had five or more individuals (Table 1). Overlap between data group UDs was calculated using Bhattacharyya's affinity 119 . The null hypothesis of no spatial difference in range use was tested by permuting the data group labels (year or stage) 1000 times and calculating overlap for each permutation (e.g. stage or year, as appropriate). The p-values for the permutation tests were estimated as the proportion of times the observed overlap was greater than the permuted overlap. Following the same procedure, the influence of down-sampling the data to one location per day was tested to quantify overlap between UDs calculated with the original (24 locations per day) and down-sampled data (one location per day) between 2017 and 2019 for premoult penguins from Dassen Island. No significant differences were found (54% UD: proportional overlap = 0.89; p-value = 1.00) and 90% UD: proportional overlap = 0.94; p-value = 1.00).
To visualise purse-seine fishing pressure within the inshore regions of South Africa's Exclusive Economic Zone, sardine and anchovy catches caught between September and December were aggregated into a grid with a 0.16° resolution for each year of the study (2012-2019, except for 2016 due to incomplete data). An average was then taken across years to represent the relative purse-seine fishing pressure. Locations, dates and tonnage of anchovy and sardine hauls (inclusive of targeted catch and bycatch) from 2012 to 2019 were taken from vessel logbooks (Department of Forestry, Fisheries and the Environment unpub. data).

Identification of marine Important Bird Areas. Following methods developed by Birdlife
International 58,59 , IBAs for the African penguin during their pre-and post-moult stages were identified using the R package track2KBA 120 . An IBA is defined as a site that is known to regularly hold significant numbers of globally threatened species or a site that supports > 1% of the global population of a congregatory bird 59,121 . IBAs were identified for each colony and stage (data groups) but only if a data group consisted of > 10 individuals (Table 1). For each data group, areas where > 10% of individuals' core areas overlapped were identified. These areas were then assessed to check if they were representative of the tracked population and thus adequately described the at-sea distribution of the data group by calculating the 'representativeness' of the data group (function: repAssess, R package track2KBA 120 ). They were deemed representative if their representativeness score was > 70%. A > 70% 'representativeness' score of the data allows for the assumption that the sampled tracks were able to adequately identify commonly used or important areas of the population. It does not however account for the different www.nature.com/scientificreports/ duration of the tracks and associated missing fixes that may result from random loss of transmission due to factors like depleted battery or device malfunction.
To enhance the practicability of management zones, spatial polygons were aggregated to minimise the boundary-to-area ratio, following the methods of Handley et al. 122 . Specifically, any isolated polygon or hole within a larger polygon, smaller than 5% of the total area identified, was removed or filled, respectively, using the R package smoothr 123 . Polygons were further merged if the great circle distance between their centroids was < 5% of the distance between the two most distant polygon centroids. The final boundaries of sites identified for each data group were delimited by a minimum convex polygon.
The relative overlap of IBAs with existing no-take zones for purse-seine fisheries 124 was calculated. Following this, the relative overlap of IBAs with the proposed 20 km no-take zones to purse-seine fisheries around Dassen, Robben, Bird and St Croix islands was calculated. To investigate the importance of these areas to pre-and post-moulting African penguins the percentage of time spent within 20 km (representative of the proposed notake zones) and 50 km (representative of the maximum foraging range of breeding penguins) 41,56 of a penguin's colony was calculated (using the percentage of fixes). Percentages were calculated both excluding and including locational fixes recorded at the colony.
Values are given as means ± standard deviations, unless otherwise specified. Significance is set at p ≤ 0.05. All data analyses were performed in the R statistical environment R version 4.0.5. 125 .

Data availability
Requests for tracking data used in this study may be made via the BirdLife International Seabird Tracking Database (http:// www. seabi rdtra cking. org).